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Non-thermal escape of molecular hydrogen from Mars 

M. Gacesa 1 ' 2 , P. Zhang 1 , and V. Kharchenko 1 ' 2 



We present a detailed theoretical analysis of a non- 
thermal escape of molecular hydrogen from Mars induced by 
collisions with hot atomic oxygen from martian corona. To 
accurately describe the energy transfer in O + H.2(v, j) col- 
lisions, we performed extensive quantum-mechanical calcu- 
lations of state-to-state elastic, inelastic, and reactive cross 
sections. The escape flux of H2 molecules was evaluated 
using a simplified ID column model of the martian atmo- 
sphere with realistic densities of atmospheric gases and hot 
oxygen production rates for the low solar activity condi- 
tions. An average density of the non-thermal escape flux 
of H2 of 1.9 x 10 5 cm -2 s -1 was obtained considering en- 
ergetic O atoms produced in dissociative recombinations of 
O2 ions. Predicted rovibrational distribution of the escap- 
ing H2 was found to contain a significant fraction of higher 
rotational states. While the non-thermal escape rate was 
found to be lower than Jeans flux for H2 molecules, the non- 
thermal escape rates of HD and D2 are significantly higher 
than their respective Jeans rates. The accurate values of 
non-thermal escape fluxes of different molecular isotopes of 
H2 may be important in analyses of evolution of the martian 
atmosphere. The described molecular ejection mechanism is 
general and expected to contribute to atmospheric escape of 
H2 and other light molecules from planets, satellites, and 
exoplanetary bodies. 

1. Introduction 

The interaction of the martian atmosphere with the solar 
radiation and interplanetary plasma results in its evapora- 
tion due to thermal (Jeans) escape and a number of non- 
thermal mechanisms. The absence of the intrinsic mag- 
netic field on Mars and its low gravitational potential make 
the martian atmosphere particularly susceptible to erosion 
[Acuna et al., 1998]. The current low atmospheric pressure 
can be explained well by extrapolating the escape rates to a 
geological time frame while accounting for a change of the 
solar activity in the past [Chassefiere and Leblanc, 2004]. 
The present day degradation of the martian atmosphere oc- 
curs mainly via non-thermal escape processes induced by ion 
charge-exhange, sputtering and ionospheric outflows driven 
by solar wind [Chassefiere and Leblanc, 2004; Johnson et al, 
2008]. The dissociative recombination (DR) of Oj is a ma- 
jor source of hot O atoms in the upper atmosphere of Mars, 
responsible for the escape of oxygen and formation of mar- 
tian hot corona [Ip, 1988; Fox, 1993; Krest'yanikova and 
Shematovich, 2005]. 
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Nascent hot O atoms collide with thermal constituents 
of the martian atmosphere and eject them from the plane- 
tary gravitational field, if a sufficient kinetic energy transfer 
occurs. Suprathermal neutral oxygen was shown to be im- 
portant for analyses of Mars' corona and the non-thermal 
escape of neutral atoms [Krest'yanikova and Shematovich, 
2005; Krestyanikova and Shematovich, 2006; Fox and Hac, 
2009]. Recent calculations of the non-thermal He escape 
from Mars carried out with accurate energy transfer param- 
eters predicted a significant He escape flux induced by hot 
O atoms [Bovino et al., 2011]. 

In this Letter we explore collisional ejection of molecules 
from the martian atmosphere. Specifically, we report the 
results of a quantum-mechanical study of the energy trans- 
fer from the hot O to H2 molecules and their subsequent 
escape. Significant computational difficulties arise from the 
fact that molecular internal rotational and vibrational de- 
grees of freedom can be excited in collisions. In addition, the 
reactive pathway leading to the production of OH molecules 
is energetically permitted. To account for the increased com- 
plexity, we have computed the cross sections for 0( 3 P) + H2 
reactive collision using fully quantum-mechanical approach. 
Kinetic theory was used to calculate the rate of energy trans- 
fer, as well as distributions of excited rotational and vibra- 
tional (RV) states of the recoiled H2 molecules. The total 
escape flux of H2 from Mars and RV distributions of escap- 
ing molecules have been evaluated for the low solar activity 
conditions. Also, we have estimated the non-thermal escape 
fluxes of HD and D2 and compared them to the correspond- 
ing Jeans escape rates. Finally, the dependence of the molec- 
ular ejection fluxes on the gravitational escape threshold is 
analyzed for conditions present on other planets, satellites, 
and exoplanets. 

2. Cross Sections and Energy Transfer 

The DR of Oj with electrons proceeds via five possible 
dissociation pathways, producing 0( 3 P), 0( 1 D), and 0( 1 S) 
[Guberman, 1988; Fox and Hac, 2009]. Energetic metastable 
0( 1 D) atoms decay via spontaneous emission and quench- 
ing in collisions with atmospheric gases into 0( 3 P) atoms 
[Kharchenko et al., 2005]. The cross sections for 0( 3 P) 
and 0( D) colliding with He were found to be very similar 
[Bovino et al., 2011]. For simplicity, we assumed a similar 
behavior for 0( 3 P) + H2 and 0( X D) + H2 elastic collisions. 

To describe the collision ejection of H2 molecules, we have 
calculated elastic and inelastic cross sections for the center- 
of-mass (CM) collision energies from 0.01 to 4.5 eV. The 
quantum scattering code ABC [Skouteris et al., 2000], that 
can treat elastic, inelastic, and open reactive channels of 
the OH production, 0( 3 P) + H 2 (v, j) -> 0~R(v",j") + H, 
where (v,j) and (v",j") indicate initial and final RV lev- 
els of H2 and OH, was used to solve the time-independent 
coupled-channel Schrodinger equation in Delves hyperspher- 
ical coordinates. In addition, for high collision energies, the 
elastic and inelastic cross sections were calculated using the 
MOLSCAT [Hutson, 1994] code to ensure the convergence of 
the nonreactive channels 1 . Extensive numerical convergence 
tests were carried out for the both codes. 
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Figure 1. Elastic and momentum transfer cross sections 
for H2(u = 0,j = 0, 1, 2) + O collisions. The momentum 
transfer cross section shown is thermally averaged over 
the first three rotational states. Inset: Total inelastic 
cross sections o™ l (E) = ^2 V , ■, a v j iV >jt (E) for H2(u = 
0, j = 0, 1,2) + O H 2 (v' = 0, j") + O, and j / j' . The 
reactive cross section for H2(v = 0, j = 0) + O — > OH + H 
is also shown and compared with the experimental results 
(black squares) [Garton, 2003]. 



The 0( 3 P)+H2(u, j) interaction was described using two 
lowest potential energy surfaces, Rogers' LEPS 3 A" [Rogers 
et al, 2000] and Brandao's BMS1 3 A' [Brandad et al, 2004]. 
Partial cross sections for initial and final rotational levels 
j and j' were constructed as a statistically weighted sum 
of the independently calculated cross sections for the two 
potential surfaces, where both 3 A" and 3 A' contribute a 
weight factor of 1/3 [Balakrishnan, 2004]. Elastic, inelastic, 
and momentum transfer partial cross sections for oxygen 
colliding with the hydrogen molecule in three energetically 
lowest rotational states are given in Figure 1. We compared 
our reactive cross sections for OH production to the previ- 
ously published results [Balakrishnan, 2003, 2004; Braun- 
stein et al., 2004; Wei-keh Wu, 2010] and found them to be 
in close agreement within the available energy range. 

To determine the energy transfer rate from the suprather- 
mal oxygen to atmospheric H2 and find its escape rate, we 
used kinetic theory with a quantum description of internal 
molecular structure and realistic anisotropic cross sections. 
Since the reactive cross sections are an order of magnitude 
smaller than the elastic cross sections (Figure 1), and the 
more massive OH molecule has a considerably higher escape 
threshold than H2, we neglected it in this study. However, 
note that a small fraction of produced OH molecules may 
be sufficiently energetic to escape. The transferred energy 
from the energetic projectile O to the frozen target H2 in the 
laboratory frame (LF) can be expressed as [Johnson, 1982] 
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where mo and mn 2 are masses of O and H2, respectively, E 



is the collision energy in the LF, j v i 
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of e v /ji and e, the CM translational kinetic energies after 
and before the collision, respectively, and 8 is the scattering 
angle in the CM frame. The energies e v iy were calculated 
quantum mechanically for the two triplet potential surfaces. 
Eq. (1) takes into account that the energy transferred to H2 
molecules is spent on increasing their translational kinetic 
energy and exciting their internal RV degrees of freedom. 



The fraction of energized H2 molecules capable of escap- 
ing can be calculated as 
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where Q v jyji (9) is the differential cross section for scat- 
tering of H2 in the initial (v,j) into the final (v',j r ) state. 
The critical angle f5 m i n was determined from the condition 
that the translational part of the transferred energy T v i i3 -/ 
is equal to the minimum energy required for H2 to escape 
from Mars, E eac = 0.26 eV. An alternative description of the 
escape process could be constructed by performing Monte 
Carlo simulations with accurate quantum cross sections for 
angular distributions of the recoiled H2 molecules. 
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collisions were calculated using [Parker and Pack, 1978] 
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3. Flux and Distribution of Escaping H 2 

Jeans escape and collisions with hot oxygen are the ma- 
jor mechanisms that contribute to the escape of neutral H2 
molecules and their isotopes from the martian atmosphere. 
Both processes are strongly dependent on the temperature 
and density of upper layers of the martian atmosphere. The 
temperature of the exosphere (above the altitude of about 
180 km), Toxo, is approximately constant and estimated to 
be between 240 and 280 K, depending on the solar activ- 
ity and gas density profiles [Krasnopolsky , 2010; Fox and 
Had, 2009; Fox, 2003]. To obtain a conservative estimate 
of the non-thermal flux of escaping H2, we considered the 
Tcxo = 240 K, corresponding to the low solar activity. Fur- 
thermore, we assumed a thermal distribution of the initial 
rotational states of H2, where more than 95 % of the to- 
tal population is distributed between its first three rota- 
tional states, j — 0, 1, 2, with corresponding population frac- 
tions equal to 0.31, 0.46, and 0.19, respectively. Using Eqs. 
(2,3) we have calculated the values of the thermally aver- 
aged momentum transfer cross sections and the fractions 
Y v v] 3 (E) of rovibrationally excited H2, sufficiently energetic 
to escape from Mars (Figure 2). The fraction Y v v ? (E) be- 
comes significant at collision energies greater than 0.7 eV 
for H2(j' = 0, 2). Although higher rotational states require 
increasingly larger projectile energies, e.g. H2(j' = 16) can 
escape for E > 1.65 eV, their fraction in the RV distribu- 
tion of the escaping molecules also becomes larger. Note 
that the initial population of higher vibrational levels of H2 
at the exobase is negligible. 

Since Y^? (E) depends only on the energy transfer ef- 
ficiency and the escape energy threshold, it can be easily 
generalized to different astronomical objects. We illustrate 
this for two hypothetic planets, the first corresponding in 
size and mass to Earth and the second to the extrasolar 
super-earth Kepler-lOb (3.3 Earth masses) [Batalha et al., 
2011]. Note that only very energetic H2, mostly in higher 
excited rotational levels, is able to escape (Figure 2). 
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Figure 2. Top: State-to-state momentum transfer cross 
sections for the H%(v' Curves for the final rotational 
levels, f = 0, 2, ... , 16 (labeled on graph), collisionally 
excited from the ground state of Eh, are shown. Bottom: 
Fraction F^ J of the recoiled H 2 (v' = 0, j') with energies 
greater than the escape energy for Mars. The escaping 
fractions for Earth for j' = 0,4 (dashed blue) and Kepler- 
10b super-earth for j' = 0, 8 (light green) are given for 
comparison. 



A simple estimate of the total escape flux of H 2 (i/,j') 
can be obtained from the exobase approximation [Fox, 2003; 
Chassefiere and Leblanc, 2004; Krasnopolsky , 2010], using 
the density of exospheric H2, and fractions T"? (E) calcu- 
lated above. However, such an approach neglects the hot O 
production below the exobase and loss of the upward flux 
in atmospheric collisions, resulting in a large uncertainty in 
the computed flux. 

We constructed a more realistic ID model of escape, anal- 
ogous to the one used to describe the escape of He atoms 
[Bovino et al, 2011]. In our model the explicit consideration 
of energy transfer collisions is combined with the altitude- 
dependent rate of production of hot O atoms, f(E, h) via DR 
channels [Guberman, 1988; Petrignani et al., 2005; Bovino 
et al., 2011]. In addition, we estimated the extinction of 
fluxes of suprathermal O and H2 due to collisions with ther- 
mal atmospheric gases. All calculations were performed for 
low solar activity. We used the rate of production of hot O 
below 400 km by Fox and Hac [2009] and smoothly inter- 
polated it to the rate given by Krasnopolsky and Gladstone 
[1996] at higher altitudes. 

The volume production rate of escaping hot H2(u', j') can 
be expressed as 



Pvy(h) = - J dETn 2 (h,E)nu 2 (h)r v v / (E)<, (j ,(£) 
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12 3' 
Volume production rate of escaping H 2 (1 3 cm 3 s 1 ) 

Figure 3. Altitude profile of the volume production rate 
Pv'.j'(h) of the non-thermal flux of H2 molecules escap- 
ing from Mars. The most significant rates with respect 
to the initial and final rotational states, j — (solid), 
j = 1 (dashed) and j = 2 (dotted), and j' = — 10, are 
shown. The curves are denoted as jj' . 



The transparency factor Tk 2 is equal to the escape prob- 
ability of hot H2(i/, j') produced in collisions with the inci- 
dent hot O of energy E at the altitude h 2 . The second trans- 
parency factor To is defined as the probability that the hot 
O atoms, produced at the altitude h 2 , reach the altitude h 
without the energy loss in collisions with other atmospheric 
constituents. The quantity n Ha (h) T v v 'f (E) o-^f yj , (E) is the 
inverse mean free path for 0+}i 2 (v, j) collisions, resulting 
in the energy transfer greater than the H2 escape threshold. 
The prefactor 1/2 indicates that, in our simplified ID model, 
approximately half of the nascent energetic atoms and re- 
coiled H2 molecules are scattered towards the planet and 
cannot escape regardless of the energy transferred. Summa- 
tions of the flux loss of H2 and O in collisions with the i-th 
atmospheric gas of density rii(h) and momentum transfer 
cross section a^ i (E) and ao^i(E), respectively, included 
major constituents of the martian upper atmosphere: CO2, 
CO, N2, O2, H2, H, Ar, and He. The momentum transfer 
cross sections for H-H2 [Krstic and Schultz, 1999], Ar-H2 
[Uudus et al, 2005], H2-H2 [Phelps, 1990] were used from 
the literature. Since no data were available in the required 
energy range, we used approximate mass-scaled cross sec- 
tions for He-H 2 (from Ar-H 2 ), N 2 -H 2 , 2 -H 2 , CO-H 2 (from 
0-H 2 ), and C0 2 -H 2 (from 0-N 2 [Balakrishnan et al., 1998]). 

Using Eq. (4) we have calculated volume production 
rate of the escaping H2(«' = 0,j') molecules induced in 
H 2 (i> = 0, j — — 2) + O collisions for a range of alti- 
tudes from ftmin = 130 km to h max = 800 km (Figure 3). 
Note that, by symmetry arguments, for a homonuclear H2 
only Aj — 0, 2, 4... transitions are allowed [Cohen- Tannoudji 
et al, 1986]. The resulting escape rates of H2 molecules are 
the largest for the elastic collisions, followed by the three 
times smaller rates for the first two excited rotational states, 
j' = 2 and j' = 4. The rates remain significant for the final 
rotational levels up to j' = 10. 
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Figure 4. Top: Collisional and thermal total escape 
rate of H2(u' = 0,j') for the first 16 rotational states j'. 
Bottom: The same as above for H 2 , HD, and D2. 



The altitude profile of the production rate of H2 capable 
of escaping is similar to the production rate profile of He 
[Bovino et ai, 2011]. This was expected, since the escape 
of both species is driven by collisions with the nascent fast 
O atoms, produced mostly below 150 km for the considered 
atmospheric and solar conditions. The calculated altitude 
profile can be used to compute the non-thermal escape flux 
of H 2 molecules and estimate the accuracy of the exobase 
approximation. We calculated </y , the non-thermal flux for 
Ha(v' = 0, j') molecules as 



,(h)dh. 



(6) 



Total collisional and thermal fluxes were calculated as sums 
over all rotational levels and found to be 1.9 x 10 s cm -2 s _1 
and 1.1 x 10 6 cm -2 s _1 , respectively. A comparison of Jeans 
and non-thermal rates of escape of H 2 (V = 0,j') molecules, 
escaping in different rotational states j' from martian day- 
side, is given in Figure 4. To simplify the calculation we 
assumed the average solar conditions and neglected the lat- 
itude dependence of the production rate of hot O atoms. 
Jeans rate is about eight times greater than the non-thermal 
rate of the escaping H 2 for the lowest three rotational states, 
while for j' > 3 the latter starts to dominate. The distinct 
character of the two RV distributions is a clear signature of 
different physical escape mechanisms. 

While, in case of H 2 , thermal rate is almost an order of 
magnitude higher than the collisionally-induced rate of es- 
cape, relative importance of the two processes changes for 
heavier isotopologues, namely HD and D 2 (Table 1). A sim- 
ilar scaling of the collisional and thermal escape fractions 
can be expected for molecular escape from more massive 
astronomical objects. 



4. Conclusions 

We find that the collisionally-induced outflow of H2 
molecules and their heavier isotopes contributes to the evo- 
lution of the martian atmosphere. Namely, the escape rate 
of molecular hydrogen induced by collisions with hot oxy- 
gen from the martian atmosphere was calculated and found 
to be about six times smaller than the corresponding Jeans 
escape rate for the low solar activity. For heavier molecules, 
the collisional escape will dominate over thermal, as we have 
illustrated in case of HD and D 2 isotope molecules. In fact, 
the described process of molecular ejection induced by colli- 
sions may be one of the most important escape mechanisms 
of HD and D2 from Mars. Consequently, the calculated es- 
cape fluxes of H 2 isotopologues could be important in anal- 
yses of the H/D ratio on Mars and evolution of water in 
martian history. 

The described mechanism of molecular escape, where col- 
lisions provides sufficient translational energy to exceed the 
escape threshold and simultaneously excite internal molecu- 
lar degrees of freedom, is rather general. It could be used to 
evaluate non-thermal escape fluxes and RV distributions of 
heavier molecules, such as CO, N 2 , or CH4, from Mars, Solar 
system bodies, and exoplanets. For the escape flux induced 
by O atoms produced in DR the upper limit on the mass of 
the escaping molecule is about 30 u on Mars. Similarly, we 
estimate that the non-thermal escape of H 2 from a planetary 
atmosphere is possible for planetary masses up to about 3.4 
Earth masses. A number of solar system bodies as well as 
the lightest currently confirmed exoplanets belong in that 
mass range. These limits do not include other non-thermal 
sources of hot atoms. 

The escaping H 2 molecules exhibit a characteristic inter- 
nal energy distribution, with a significant fraction of pop- 
ulated higher rotational states. Since H 2 molecules do not 
have a permanent electric dipole moment, they decay to the 
ground state mainly via collisions with the background gases 
present in an extended planetary corona. This is true for all 
escaping molecules that do not have ther permanent dipole 
moment. It could be possible to indirectly detect the pres- 
ence of H 2 or other rotationally excited light molecules in the 
extended martian corona from a careful analysis of the col- 
lision rates and abundancies of the excited coronal species. 
Finally, a significant amount of rovibrationally excited H2 
molecules remain in the martian atmosphere after colliding 
with hot O atoms. The cross sections and energy transfer 
parameters presented in this study can be used to determine 
non-thermal translational and rovibrational distributions of 
hot H 2 gas in the upper atmosphere of Mars. 
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Notes 

1. A detailed description of the scattering calculations and result- 
ing cross sections for the two lowest potential energy surfaces 
will be published elsewhere. 



Table 1. Total collisionally-induced escape rates of H2, HD, 
and D 2 from the martian atmosphere. 
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